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ABSTRACT 

The  flow  of  the  classical  linearly  viscous  fluid  between  two  infinite 
parallel  planes  rotating  with  constant  (but  different)  angular  velocities 


about  a  common  axis  has  received  a  great  deal  of  attention  during  the  past  60 
years  (cf.  Parter  tlzTVi  However,  until  recently  the  emphasis  has  been  on 
solutions  which  are  axi-symmetric.  Recently,  Berker  [3]  in  his  study  of  the 
flow  between  parallel  planes,  rotating  with  the  same  constant  angular 
velocities  about  a  common  axis  exhibited  a  one  parameter  family  of  solutions 
only  one  of  which  is  axi-symmetric.  In  this  study ,  -we  exhibit  the  existence 
of  a  one  parameter  family  of  solutions  (for  -large'^  viscosities)  when  the 
planes  are  rotating  with  constant  but  different  angular  velocities  about  a 


common  axis  or  about  non-coincident  axes. 


AMS  (MOS)  Subject  Classifications:  34B15,  34D10,76D05 

Key  Words:  Ordinary  differential  equation,  rotating  fluids,  non-uniqueness 
Work  Unit  Number  1  -  Applied  Analysis 
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SIGNIFICANCE  AND  EXPLANATION 


The  flow  between  parallel  planes  rotating  with  constant  (but  different) 
angular  velocities  about  a  common  axis  was  first  studied  by  Batchelor  [2] 
instigated  by  an  earlier  study  of  von  Karman  [9]  in  1923  and  it  has  been  the 
object  of  considerable  analysis  ever  since  then.  The  Karman  assumptions, 
which  has  always  been  employed  in  studying  the  problem,  leads  to  solutions 
which  are  axi-symmetric.  In  this  study  we  show  that  the  classical  solutions 
are  imbedded  in  a  one  parameter  family  of  solutions  which  are  basically  non- 
symmetric  and  hence  the  classical  solutions  are  not  "isolated"  or  "stable". 
These  results  which  seem  to  arise  due  to  the  unboundedness  of  the  flow  domain 
lead  one  to  question  the  meaning  and  validity  of  the  study  of  the  above 
boundary  value  problem.  Similar  one  parameter  family  of  solutions  are  shown 


to  exist  in  the  case  of  the  flow  of  rotating  planes  of  constant  but  differing 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
sumsiary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


REMARKS  ON  THE  FLOW  BETWEEN  TWO 


PARALLEL  ROTATING  PLATES 


Seymour  V.  Parter  and  K.  R.  Rajagopal 


Introduction 

Let  II,  and  Ilj  be  two  infinite  planes  parallel  to  the  (x,y)  plane,  say  nl  is 
the  plane  s  ■  -1  while  IIj  is  the  plane  t  m  1.  Let  a  >  0  be  a  fixed  constant  and 
suppose  that  II ^  rotates  about  a  point  (x  «  0,  y  ■  -a/2,  z  «  -1)  with  constant  angular 
velocity  (I  f  while  II2  rotates  about  the  point  (x-0,  y  -  +a/2 ,  i  »  1)  with  constant 
angular  velocity  fl+. .  We  suppose  a  classical  incompressible  fluid  fills  the  infinite 
space  between  these  planes  and  we  seek  steady  state  solutions  of  the  Navier-Stokes 
equation  which  describe  the  fluid  flow. 


It 


-1 


Finally,  we  make  the  basic  assumption  that 
1.1)  Um  -  -H(z>, 

that  is,  the  component  of  velocity  in  the  z  direction  is  a  function  of  z  alone. 

If  a  -  0  and  we  also  assume  that  the  flow  is  axi-svmmetric  then  the  basic  theory  of 
von  Karman  [9]  and  Batchelor  [2]  leads  to  the  following  conclusions: 
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a 


1.2) 

0  -§H‘(a) 

r  2 

1.3) 

°e  “  2  6(*) 

where  the 

function*  <G(z),H(z)> 

1.4a) 

1 V 

CH  +  HH' * '  +  GG* 

1.4b) 

EG"  +  HG1  -  H'G  - 

r  -  (xV  >4 

are  solutiona  of  the  boundary-value  problem 
-  0,  -1  <  a  <  1 
0,  -1  <  *  <  1 


and 

1.5a)  H(-1,e)  «  H(1,e)  -  0  (no  penetration) 

1.5b)  H'  (-1,e)  -  H'(1,e)  -  0  (no  elip) 

1.5c)  G(-1,e)  -  20_1  G ( 1 , e )  -  20+1 

where  the  poaltlve  parameter  c  >  0  la  related  to  the  bulk  viscoalty 

Thle  boundary  value  problem  haa  been  atudled  at  great  length  There  are  many 
numerical  atudlea  and  many  formal  asymptotic  studlea.  There  are  aleo  rlgoroua  exletence 
theorema 


(I)  for  e  >>  1  by  Haatlnga  [8]  and  Herat  [6], 

(II)  for  n_1  *  -fl+1  +  0  and  all  e  >  0  by  McLeod  and  Parter  (11), 

(III)  for  0  <  e  «  1  by  Krelaa  and  Parter  [10]. 

The  recent  aurvey  article  [12]  containa  a  reaaonably  up-to-date  dlacuaalon  of  thla 

problem. 

When  0_1  -  Q+1  +  Q  there  la  one  apeclal  aolutlon  (not  the  only  aolutlon-aee  [3])i 
the  rigid  body  rotation  given  by 
1.6)  HU,e)  s  0,  G(«,e)  s  2Q_^ . 

It  la  not  difficult  to  verify  ([5],  [18])  that  thla  aolutlon  la  "atable"  and  "iaolated" 
relative  to  the  von  Karmen  equatione  (1.4),  (1.5).  By  "iaolated*  we  mean  there  la  a 
neighborhood  of  thla  aolutlon  wherein  there  are  no  other  aolutiona,  and  by  "atable"  we 
mean  there  la  no  bifurcation  from  thla  aolutlon.  In  particular  the  linearixed  problem  at 
thla  aolutlon  la  non-aingular.  On  the  other  hand,  Berker  [3]  haa  conatructed  a  one 
parameter  family  of  aolutiona  to  the  general  ataady  Havier-Stokea  equation  which  includea 
the  rigid  body  motion.  The  rigid  body  motion  la  the  only  axiaymmetric  aolutlon  that 
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belongs  to  Barker's  [3]  family.  In  Cartesian  coordinates  his  solution  takes  the  fore  (we 


have  set  Q_1  *  fl+1  «  1). 

1.7)  0  »  -[y-g(z)],  0  *  [x-f(z)),  U  -  0 

x  y  z 

where 

1.8a)  f(i)  -  t  [♦(z)-*(a)l  -  tx(z>-x<1>)}, 

1.8b)  g(z)  -  t  (♦(z>-*(1)]  +  (x(z)-x(l))}. 

with 


1.9a) 

1.9b) 

1.9c) 


+(z)  -  coshmz  •  cosmz. 


X(z)  -  sinhaz  •  sinmz, 

,  i  ,  v2 

"  “  ‘ie’  2 


1.9d)  A  “  [  1  — ♦  ( 1 ) )  2  +  (xd)]2  “  (coshm-cosm)2 

while  t  is  an  arbitrary  positive  constant.  Observe  that  (1.7)  shows  that  this  solution 
satisfies  the  basic  ansatz  (1.1). 

The  case  a  f  0  and  fl_1  *  Q+1  «  0  has  relevance  to  the  flow  occurring  in  the 
orthogonal  rheometer,  an  instrument  that  is  employed  in  determining  the  material  moduli 
which  characterize  non-Newtonian  fluids.  Recently.  Rajagopal  [13]  has  studied  the  flow  of 
general  simple  fluids  in  such  a  domain  and  Rajagopal  and  Gupta  [15] ,  and  Rajagopal  and 
Wineman  [17]  have  established  exact  solutions  to  the  problem  for  certain  non-Newtonian 
fluids?  Rajagopal  and  Gupta  [18]  have  also  established  a  one  parameter  family  of  exact 
solution  for  an  incompressible  homogeneous  fluid  of  second  grade  when  a  •  0  and 


a_,  -  n+1  -  n. 

When  a  *  0  and  -  fl+1,  the  solution  corresponding  to  the  usual  Kerman  ansatz 

leads  to  exactly  one  of  the  solutions  in  Barker's  [3]  one  parameter  fairly,  namely  the 
rigid  motion  which  is  axi-symmetric.  In  the  case  a  +  0,  and  “  Q+1,  an  exact 

solution  hss  been  obtained  for  the  classical  incompressible  fluid  by  Abbot  and  Halters 


More  recently,  Goddard  [7]  has  also  established  results  which  are  similar  to  those  in 
[131. 


1TTT 
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[1].  The  existence  of  such  a  solution  Is  lulled  In  an  earlier  analysis  of  Barker  [4). 

This  motivates  us  to  look  for  a  more  general  class  of  solutions  to  the  Karman  problem  (and 

to  the  corresponding  problem  when  a  f  0)  which  would  reduce  to  the  class  of  non- 

axisymmetric  solutions  exhibited  by  Barker  [3]  Thus,  we  seek  a  solution  -  for  classical 

incompressible  fluids  -  which,  in  Cartesian  coordinates  -  take  the  form 

1.10a)  Ux  -  |  H'(s)  -  y  +  g<«) 

1.10b)  Vy  “  2  H'(*>  +  x  "  f(l) 

1.10c)  U  -  -H(s). 
s 

In  cylindrical  coordinates  this  velocity  field  takes  the  for* 

1.11a)  vr  -  f  H*  (z)  +  g(s)  coe  0  -  f(s)  sin  6, 

1.11b)  v#  “  f  G<*>  "  9<*>  *in  e  "  f(B)  “•  e> 

1.11c)  V  -  -H(s)  . 

z 

Observe  that,  if  H  -  0  and  G  -  2fl  we  have  a  velocity  field  of  the  form  described  by 

Barker,  And,  if  f  -  g  ”  0  we  have  a  velocity  field  of  the  form  described  by  von  Karman. 

As  we  shew  in  section  2,  there  is  a  solution  of  the  steady  state  Navier-Stokes 
equations  of  the  form  (1.10)  (or  (1.11)]  if  and  only  if  the  functions  H(z),  G(z),  g(z), 
f(s)  are  solutions  of  the  boundary-value  problem 


1.12a) 

Av 

EH  +  HH*  ' '  +  GG'  -  0 

,  -1  *  •  *  1 

1.12b) 

CG"  +  HG'  -  H'G  -  0, 

-1  <  z  «  1 

1.13a) 

ef'"  +  H If"  ♦ViH'f' 

♦!$  (Gg)'  -  0 

1.13b) 

eg"'  +  Hg"  ♦  VjjH'g' 

-Vi""  9  -'t  («)’  -  o 

1.14a) 

■  (-1,0  -  8(1, t)  -  0, 

(no  penetration). 

1.14b) 

■  •(-1,0  -  H'(1,e)  -  0 

,  (no  slip). 

1.14c) 

G(-1,c)  -  2Q_1 ,  G(1,e) 

-2fl+1 

1.15) 

f(-1,e)  -  f(1,e)  -  0, 

g(-1,e)  -  -afl^/2,  9(1, e)  -  a(l+1/2 

Me  note  that  equations  (1.12)  with  boundary  conditions  (1.14)  are  exactly  the 
nonlinear  von  Karman  equations  for  axially  symmetric  swirling  flow  for  functions 
<  H(s,e),  Q(s,c)  >  while  the  equations  (1.13)  with  boundary  conditions  (1.15)  are  linear 


aquations  for  <  f(z,e),  g(z,e)  >  with  coafficients  depending  on  <  H(z,e),  G(z,c)  > 
which  reflect  the  asyenetry  and  the  poaaible  displacement  of  the  centers  of  rotation  of 
the  bounding  planes.  Moreover,  given  <H(z,  e) ,  G(z,e)>  the 

of  two  third  order  equations  with  only  four  boundary  conditions.  Whenever,  there  is  a 
solution  of  the  von  Karman  equations  one  can  ask  two  questions: 

(i)  In  the  case  when  a  ■  0,  is  the  axi-aymmetric  flow  inbedded  in  a  continuous  one 
parameter  family  of  sure  general  solutions? 

(ii)  When  a  M,  does  this  axi-symnetric  flow  fora  the  basis  for  a  one  paraaeter 
faaily  of  solutions  of  the  problem  for  rotations  about  different  centers? 

In  case  (i)  this  is  a  homogeneous  underdetermined  system  and  the  answer  is  veal  We  need 
merely  consider  the  additional  condition 

1.16)  f • (-1 )  -  g' (-1 )  -  0. 

If  this  augmented  homogeneous  problem  has  a  non-trivial  solution 
<f(c,c),  g(x,e)>  +  <0,0>,  then  <if(z,c),  tg(z,e)>  is  also  a  solution  for  every  real 
number  t.  On  the  other  hand,  if  the  system  (1.13),  (1.15),  (1.16)  does  not  have  a 
nontrivial  solution;  then  the  problem  given  by  (1.13),  (1.15)  and 

1.17)  f*(-1)  -  i,  g*(-1)  -  0 

yields  a  unique  solution  <f (z, e, 1) ,  g(z,c,t)>  of  the  form 

<f(z,c,l),  g(z,e,t)>  •  <tf(z,e,1),  ig(z,e,1)>. 

This  simple  result  has  the  following  important  consequence.  In  the  classical  case  of  two 
infinite  parallel  planes  rotating  about  a  cowon  axis,  (i.e.,  a  -  0)  whenever  there  is  a 
solution  of  the  von  Karman  equations  (12),  (14),  this  axi-symmetric  flow  is  inbedded  in  a 
one-parameter  family  of  solutions  of  the  full  Navier-Stokes  equations.  Thus,  despite  the 
intense  interest  in  the  von  Karman  problem  within  the  class  of  all  solutions  of  the 
Havier-Stokea  equations,  these  special  solutions  ars  "unstable”.  While  it  seems  likely 
that  a  similar  simple  argument  settles  the  natter  for  case  (ii),  it  is  not  apparent. 

In  either  case  one  can  ask  a  more  subtle  question:  can  we  find  a  family 


<f(z,e,t),  g(z,c,i)>  which  is  continuous  in  both  e  and  t  and  (at  the  same  tine)  has 


the  geometric  significance  of  the  Barker  solution  for  the  special  case  8-1  "  ♦  0? 

In  other  words,  can  we  find  solutions  of  (1.12),  (1.13),  (1.14),  (1.15)  and 
1.18)  g(0,e,t)  -  0,  f(0,e,t)  -  t. 

Since  the  system  (1.13),  (1.15)  is  linear,  the  answer  is  'yes'  for  problem  (i)  if  and 
only  if  it  is  also  'yes'  for  problem  (ii). 

In  section  3  we  answer  these  questions  in  the  affirmative  for  large  e.  While  this 
result  is  an  immediate  consequence  of  the  implicit  function  theorem  (applied  at  c  “  “) 
we  will  give  a  complete  proof. 


2.  equations  of  Motion 

In  this  discussion  we  follow  the  outline  of  the  argument  given  in  [12,  section  2].  A 
velocity  field  of  the  form  (1.10)  satisfies 

2.1)  div  U  -V^HMs)  +1/iH'd)  -  8'(t)  -  0. 

Thus,  the  basic  constraint  of  the  Havier-Stokes  equations  is  satisfied.  We  new  turn  to 
the  aquation 

2.2)  pAu  -  pu* Vu  -  Vp 

where  u  denotes  the  viscosity,  p  the  density,  and  p  is  the  pressure.  We  eliminate 
the  pressure  by  taking  the  curl  of  both  sides  of  (2.2)  and  obtain 

2.3)  liAw  -  p(umD)  “  0 

where 


2.3b) 


u 


curl  0. 


A  detailed  calculation  now  yields 

I  •  •  £1  I  I 


2.4) 

while 


Au 


-i{j  H 


1V  ♦  f  G 


>  +  1  H  -  ?  G  +  f  >  +  k  {G-} 


*  V  (6a)'  (fn')'  (f'H)' 

«  X  0  ■  -i  {|  l(GH')'-  (G'H) '  ]  -  |  [GG'  +  HH">)  +  +  Y  )  ) 


-6- 


(Gf  >'  (<*■)'  (q'H)' 


2-5)  +  i  {  -  |  IGG*  +  HH"  ']  +  |  (tG^)*  -  (GH')M  +•  ^ 


1  > 


+  k  {  GH'  -  HG'} 


On  equating  the  coefficients  of  k_  in  (2.3a)  we  obtain 

2.6a)  11  G"  +  HG'  -  H'G  -  0. 

P 

On  equating  the  coefficients  of  _i  x  we  obtain 

M  G'"  +  (HG')'  -  (H'G)'  -  0 
P 

which  also  follows  from  (2  6a).  On  equating  the  coefficients  of  i  y  we  obtain 

2.6b)  ^  HiV  +  HH*  1  1  ♦  GG'  -  0. 

P 

Finally,  the  zeroth  order  terns  in  the  coefficient  of  i.  yield 
2.6c)  p  f'"  +  "!$  (H'f  )  *  +Vfc  (Gg)’  -  0. 

The  coefficients  of  Jx,  jy  and  J_  yield  the  final  equation 
2.6d)  p9‘"  +  «9')'  “^(H'g)’  +1/i<Gf)1  -0. 

Thus  we  have  established  the  equations  (1.12)  and  (1.13)  with 


We  now  turn  to  the  boundary  conditions.  For  our  isswdiate  purposes  it  would  be 
sufficient  to  lepoee  the  conditions  (1.4)  and  (1.5)  and  observe  that  the  velocity  field 
(1.10)  now  satisfies  the  steady-state  Navier-Stokes  equations  and  the  boundary  conditions 
of  our  problan.  However,  it  is  somewhat  wore  satisfying  to  proceed  as  follows.  Equation 


xhua,  we  have  obtained  the  boundary  conditiona  (1.14).  Now,  using  (1.11b)  and  letting 
r  ♦  0  with  a  judicioua  choice  of  9  we  obtain  the  boundary  conditions  (1.15). 


3.  Existence  for  large  e  >>1 

In  this  section  we  present  what  is  essentially  a  standard  argument  for  regular 
perturbation  problems.  The  argument  is  given  in  some  detail  because  we  wish  to  emphasize 
the  following  facts. 

(i)  There  is  an  e  >>  1  and,  for  all  e  >  e  there  is  a  solution  of  the 

o  o 

von  Kansan  problem  (1.4),  (1.5).  Moreover,  this  solution  is  continuous  in  e.  Hence, 

th'jre  is  a  curve  of  solutions  and  there  it  no  local  bifurcation  of  solutions  of  the  von 

Xarman  equation  from  this  curve.  Again,  within  the  set  of  solutions  of  the  .  » \rman 

equations,  for  fixed  e  >  e  ,  each  of  these  solutions  is  isolated. 

o 

(ii)  Nevertheless,  in  this  same  range  of  e  there  is  a  solution  of  full  system 

(1.12)-(1.15):  a  one  parameter  family  of  solutions  <H(x,c),  G(x,e),  f(x,e  )  x,e,t)> 

which  includes  (for  £  *  0)  the  axi -symmetric  von  Karman  solution.  Moreover,  if 

e  and  £  are  both  fixed,  this  solution  is  an  isolated  stable  solution.  Of  course, 
with  a  >  0  these  solutions  provide  a  one  parameter  family  of  solutions  of  the  problem  of 
rotation  about  different  centers. 

Our  first  goal  is  to  show  that  in  the  case  of  the  von  Karman  equations  a  relatively 
simple  Picard  iteration  scheme  converges  for  e  >>1  and  -  in  the  nature  of  things,  the 
eolutiona  so  obtained  are  continuous  in  R  -  1/e. 

Definition!  Let  f  «  C*  1-1,1)  ,  k  >  0.  Then 

Ifl  -  l  Max{|f(j)(x)|:  -Itxtl) 
j-0 

Lemma  3. 1 :  Consider  the  two  boundary  value  problems: 

3.1a)  Ql4)  -  f,  -1  4  x  4  1 
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3.1b) 


euu  -  e*(*i)  -  o 


3.2a)  +’  -  g.  -1  <  x  <  1  , 

3.2b)  *<-1)  -  20_1  ,  *<1)  -  20+). 


There  is  a  constant  K  >1  such  that 
o 

3.3a)  I0*4  <  VflQ 

3.3b)  ***2  <  V,9,o  +  2|Q-l'  +2|fl+in 

Proof:  Direct  integration. 


Let 

3.4a) 

3.4b) 


a_t  and  fl+1  be  given.  Let 

a  -  20Ko(|a_1l+|n+1ll  +  1  , 

R  "  e  *  Ro  “  ^6K~o^ 
o 


Let 

3.5a)  H  3  0, 

o 

3.5b)  Go  -  2  ta..,  ♦  ^7^- 
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which  la  continuous  in  e  and  i.  Thus  wa  have  verified  all  tha  opening  remarks  of  this 
•action. 

4.  Remarks 

Wa  conclude  this  analysia  by  asking  a  few  observations  on  the  significance  of  the 
result  established  in  the  previous  ssetion. 

We  have  studied  special  solutions  of  the  Naviar-Stokes  equations  for  a  fluid 
contained  with  two  infinite  parallel  planes  each  rotating  with  a  constant  angular  velocity 
IJ^Ck  -  +1)  .  The  centers  of  rotation  may  or  aay  not  lie  on  the  sane  axis  perpendicular 
to  the  planes. 

In  either  case  wa  are  led  to  a  systen  of  ordinary  differential  equations  which 
contain  (as  a  subset)  tha  classical  equations  of  von  Karman  [9]  and  Batchelor  [2]  for 
special  axi-aymawtric  flow  about  a  cosoon  axis.  In  particular,  in  this  classical  case 
studied  by  von  Karman  and  Bathelor,  if  there  are  such  special  solutions,  they  are  never 
Isolated  solutions  when  considered  with  the  scops  of  the  full  Nevier-Stokes  equations.  In 
the  case  of  'off-centered*  rotation  there  are  many  unanswered  questions.  However,  we  have 
shown  that  (contrary  to  most  intuitive  ideas)  in  the  case  of  'large*  viscosity,  there  are 
solutions  and  they  are  never  isolated.  While  the  underlying  basis  for  these  anomalies  is 
not  completely  understood,  we  believe  it  is  related  to  the  fact  that  in  this  unbounded 
domain  the  velocities  at  large  r  are  great. 

It  is  also  worth  observing  that  similar  results  can  be  established  in  the  case  of  the 
flow  of  a  Newtonian  fluid  between  rotating  porous  disks  and  e  >  >  1.  In  this  case,  the 

only  change  in  the  problem  is  in  the  boundary  condition  (1.14a).  It  is  an  easy  matter  to 
modify  the  arguments  of  Section  3. 
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